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6.1  Introduction. 


6.1.1  Nonlinearity  in  Periodic  Phononic  Structures: 

The  effect  of  structural  periodicity  on  wave  propagation  has  been  studied  in  a  wide  array 
of  fields.  This  includes  electrons  in  crystalline  lattices,  light  waves  in  photonic  periodic 
structures,  cold  atoms  in  optical  lattices,  and  plasmons  in  networks  of  Josephson 
junctions  or  metal  surfaces  [1-3].  The  preceding  chapters  of  this  book  have  considered  in 
particular,  the  effect  of  structural  discreteness  and  periodicity  on  the  propagation  of 
phonons,  sound,  and  other  mechanical  waves.  Phononic  crystals  and  acoustic 
metamaterials  are  examples  of  materials  designed  for  this  purpose.  Studying  the  linear 
response  of  these  systems,  many  common  properties  are  revealed,  such  as  the  existence 
of  band  gaps,  which  influence  the  transport  properties.  However,  as  the  amplitude  of  the 
wave  excitation  is  increased,  the  response  of  the  material  becomes  nonlinear  and  the 
wave  propagation  more  complex.  The  study  of  nonlinearity  in  periodic  structures  has 
revealed  unique  phenomena  with  no  analogs  in  the  linear  theory.  Such  phenomena 
include  nonlinear  resonances,  bifurcations,  chaos,  self  trapping  and  dynamical 
localization.  Nonlinear  devices  thus  offer  potential  for  novel  applications  such  as 
frequency  conversion,  energy  harvesting,  and  switching,  among  others. 

Although  the  role  of  nonlinearity  has  been  extensively  studied  in  periodic  structures  and 
metamaterials,  such  as  photonic  periodic  structures,  optical  metamaterials  and  atomic 
Bose-Einstein  Condensates  in  optical  lattices  [4],  there  are  thus-far  few  examples  of 
nonlinear  phononic  crystals  or  nonlinear  acoustic  metamaterials  that  have  been  developed 
in  practice.  The  sources  of  nonlinearity  in  phononic/acoustic  materials  can  be  categorized 
into  (i)  intrinsic  and  (ii)  extrinsic.  The  former  derives  from  nonlinearities  in  the  material 
fundamental  constitutive  response  (i.e.,  inter-atomic  forces,  nonlinear  elasticity, 
plasticity,  or  ferroelasticity)  [5],  The  latter  derives  from  the  geometry  or  topology  of  the 
fundamental  building  blocks  (i.e.,  contact  forces  between  particles  [6],  defonnation  of 
micro-nano  mechanical  oscillators  [7],  or  the  nonlinearity  related  to  geometrical 
instabilities  [8]). 

Homogenous  materials  with  nonlinear  elastic  [5]  or  nonlinear  acoustic  response  [9]  have 
long  been  studied.  Nonlinear  bulk  and  surface  waves,  resulting  from  the  interplay 
between  the  intrinsic  nonlinearity  and  geometrical  dispersion,  have  also  been  studied  and 
observed  in  solids  [10-12],  However,  until  recently,  this  research  has  not  been  combined 
with  the  new  capabilities,  enabled  by  the  study  of  linear  phononic  crystals  and  acoustic 
metamaterials,  as  described  in  the  previous  chapters  of  this  book.  The  far  most  prominent 
example  of  nonlinear  periodic  phononic  structures  in  the  acoustic  regime  (0-20  kHz)  is 
granular  crystals,  which  are  arrays  of  elastic  particles  in  contact  [13]  whose  nonlinearity 
results  from  the  geometry  of  adjacent  particles.  In  the  ultrasonic  regime  (greater  than  1 
MHz),  nonlinear  energy  localization  has  been  observed  in  micromechanical  oscillator 


arrays  [14].  Moreover,  recent  work  by  Liang  B  et  al.  theoretically  suggested  [15]  and 
later  demonstrated  experimentally  [16]  the  ability  to  use  nonlinear  acoustic  materials, 
e.g.,  a  contrast  agent  microbubble  suspension,  coupled  to  a  linear  superlattice  to  obtain 
acoustic  rectification.  Finally,  at  much  higher  frequencies  (greater  than  1  GHz),  several 
studies  have  explored  mechanical  wave  propagation  in  periodic  nonlinear  structures, 
focusing  on  high  amplitude  stress  wave  propagation  and  the  propagation  of  thennal 
phonons.  Maris  and  collaborators  studied  the  propagation  of  high  amplitude  picosecond 
pulse  propagation  in  crystalline  solids  -  a  naturally  occurring  nonlinear  lattice  [17,  18]. 
With  respect  to  the  propagation  of  high  frequency  thermal  phonons,  many  studies  have 
focused  on  the  use  of  nonlinear  lattices  for  thennal  rectification,  the  earliest  of  which 
were  conducted  by  Terraneo  et.  al.  in  2002  [19]  and  Li  et.  al.  in  2004  [20],  An 
experimental  study  by  Chang  et.  al.,  2006,  also  demonstrated  thermal  rectification  using 
mass-loaded  carbon  and  boron-nitride  nanotubes,  and  attributed  the  rectification  to 
nonlinear  processes  [21].  Based  on  these  studies,  several  following  works  have  extended 
this  concept  further  to  suggest  nonlinear  lattices  for  use  as  transistors  [22],  logic  gates 
[23],  and  memory  [24].  Computational  studies  have  also  suggested  several  experimental 
setups  for  creating  thennal  rectification  building  blocks,  including  carbon  nano-cones 
[25]  and  graphene  ribbons  [26]. 

Because  many  nonlinear  periodic  structures  can  be  modeled  as  nonlinear  lattices,  the 
study  of  nonlinear  lattices  can  offer  many  potential  lessons  and  insights.  As  such,  the 
section  directly  following  gives  a  brief  history  of  the  major  types  of  nonlinear  lattices. 
The  review  of  nonlinear  lattices  is  then  followed  by  an  introduction  to  granular  crystals, 
which  is  one  of  the  most  widely  studied  types  of  nonlinear  periodic  phononic  structures, 
and  the  subject  matter  which  comprises  the  focus  of  this  chapter. 


6.1.2  Nonlinear  lattices 

Since  the  first  computational  experiments  in  nonlinear  mass-spring  lattices  by  E.  Fermi, 
J.  Pasta,  and  S.  Ulam  in  1955  [27],  there  has  been  a  wealth  of  interest  in  the  dynamics  of 
nonlinear  lattices.  Using  one  of  the  first  modern  computers,  Fermi,  Pasta,  and  Ulam 
(FPU)  studied  a  system  where  the  restoring  (spring)  force  between  two  adjacent  masses 
was  nonlinearly  related  to  the  relative  displacement  between  masses,  and  investigated 
how  long  would  it  take  for  long-wavelengthlength  oscillations  to  transfer  their  energy 
(thermalize)  into  an  equilibrium  distribution.  Instead  of  the  predicted  thermalization,  they 
found  that  over  the  course  of  the  simulation,  most  of  the  energy  had  returned  to  the  mode 
with  which  they  had  initialized  the  system  in  coherent  fonn. 

This  discovery  initiated  fields  of  research  relating  to  the  study  of  nonlinear  waves  in 
discrete  lattices  [28-30].  This  includes  many  different  types  of  nonlinear  lattices  inspired 


by  physical  systems  (in  addition  to  the  FPU  lattice),  and  the  study  of  physical  phenomena 
occuring  in  them.  The  nonlinear  lattices  most  commonly  studied  can  be  roughly 
categorized  into  three  types:  the  discrete  nonlinear  Schrodinger  (DNLS),  the  Klein- 
Gordon  (KG),  and  the  FPU  lattices.  The  ID  fonns  of  these  lattice  equations  are  as 
follows: 

The  DNLS  can  be  written  as 


ju,  =  -e(uM  +  km)-  |  ut  |2  ut , 

the  KG  as 

«/  =  e(w/+i  +  «w  -  2 u,)  -  V'(u,) , 

and  the  FPU  as 


ui=V'(uM-ui)-V'(ui-ui_l), 

where  w,  is  the  dynamical  variable  of  interest  at  site  i ,  e  is  a  coupling  parameter 

(constant),  j  =  V- 1 ,  and  V  is  a  nonlinear  potential  function.  The  DNLS  equation  has 
been  used  to  describe  nonlinear  waveguide  arrays  and  Bose-Einstein  condensates,  among 
others  [30].  Additionally,  under  small-amplitude  assumptions,  the  DNLS  can  be  derived 
from  the  KG  and  FPU  lattices  [31].  The  KG  system  has  been  used  to  model  systems  of 
coupled  pendula,  electrical  systems,  and  metamaterials  with  split  ring  resonators,  among 
others  [30].  In  contrast  to  the  KG  system,  the  FPU  has  no  on-site  potential  term,  and 
instead  involves  a  nonlinear  potential  based  on  nearest  neighbor  interactions  (nonlinear 
springs).  The  FPU  system  has  been  used  to  describe  the  behavior  of  crystalline  solids. 
Granular  crystal  systems  are  a  type  of  FPU  lattice. 

Studies  of  all  these  lattices  have  showed  the  emergence  of  localized  nonlinear  structures 
and  have  been  used  to  understand  the  existence  of  such  phenomena  in  other  nonlinear 
(not  necessarily  discrete)  systems.  Two  examples  of  nonlinear  coherent  structures,  which 
are  particularly  applicable  to  the  study  of  granular  crystals,  are  solitary  waves  and 
discrete  breathers.  Solitary  waves  were  first  observed  by  J.  Russel  in  a  shallow  water- 
filled  canal  in  1844  [32],  Since  then  they  were  shown  to  be  a  solution  of  the  Korteweg-de 
Vries  (KdV),  a  nonlinear  partial  differential  equation,  and  have  been  discovered  in 
myriad  systems  and  discrete  nonlinear  lattices  of  all  the  above  types  [33],  [34]  (including 
granular  crystal  systems  [13]).  Discrete  breathers  are  a  type  of  intrinsic  (not  tied  to  any 
structural  disorder)  localized  mode,  and  have  been  the  subject  of  many  theoretical  and 
experimental  investigations  [31,  35,  36],  Discrete  breathers  have  been  demonstrated  in 
charge -transfer  solids,  superconducting  Josephson  junctions,  photonic  crystals, 
biopolymers,  micro  mechanical  cantilever  arrays,  and  more  [31].  In  addition  to  nonlinear 
localized  structures,  the  presence  of  nonlinearity  in  dynamical  lattices  makes  available  an 
array  of  useful  phenomena  including  quasiperiodic  and  chaotic  states,  sub-  and 


superharmonic  generation,  bifurcations,  the  breaking  of  time-reversal  symmetry,  and 
frequency  conversion  [37-42]. 


6.1.3  Introduction  to  Granular  Crystals 

Granular  crystals,  defined  as  ordered  aggregates  of  elastic  particles  in  contact  with  each 
other,  are  a  type  of  nonlinear  periodic  phononic  structure  (figure  1).  Their  nonlinearity 
emerges  from  two  characteristics:  (1)  the  geometry  of  the  particles  is  such  that  the  stress 
transmission  at  the  contact  between  neighboring  elements  is  nonlinear;  and  (2)  in  an 
uncompressed  state,  granular  crystals  cannot  support  tensile  loads,  effectively  creating  an 
asymmetric  potential  between  neighboring  elements.  An  unusual  feature  of  granular 
crystals  that  results  from  these  nonlinearities  is  the  negligible  linear  range  of  the 
interaction  forces  between  neighboring  particles  in  the  vicinity  of  a  zero  compression 
force.  This  results  in  non-existent  linear  sound  speed  in  the  uncompressed  material,  that 
led  to  the  introduction  of  the  concept  of  “sonic  vacuum”  -  a  medium  where  the  traditional 
wave  equation  does  not  support  a  characteristic  speed  of  sound  [13]. 

The  study  of  granular  crystals  emerged  in  1983  with  the  study  by  Nesterenko,  showing 
analytically  and  numerically  [43],  and  later  experimentally  [44],  the  concept  of  “sonic 
vacuum”  and  the  fonnation  and  propagation  of  highly  nonlinear  solitary  waves  in  one¬ 
dimensional  granular  crystals.  Granular  crystals  have  since  been  shown  to  support  unique 
dynamic  phenomena,  such  as  a  tunable  response,  encompassing  linear,  weakly  nonlinear 
and  highly  nonlinear  behaviors  that  can  be  controlled  by  the  application  of  a  variable 
static  load  [13,  45-47]. 

In  their  linear  and  weakly  nonlinear  regime,  granular  crystals  have  shown  the  ability  to 
support  tunable  acoustic  band  gaps  [48,  49]  and  discrete  breathers  [50,  51].  In  the 
strongly  nonlinear  regime,  they  have  been  shown  to  support  compact  solitary  waves, 
anomalous  reflections  [52]  and  energy  trapping  phenomena  when  interacting  with  defects 
and  interfaces  [53].  The  tunability  of  their  dynamic  response  and  their  rich  dynamics 
which  has  been  confirmed  by  theory,  numerical  simulations  and  simple  experiments, 
have  made  granular  crystals  one  of  the  most  studied  examples  of  nonlinear  lattices.  They 
have  also  been  proposed  and  designed  for  use  in  numerous  engineering  applications 
including  tunable  vibration  filters  [49,  54],  optimal  shock  protectors  [55],  non-destructive 
evaluation  devices  [56],  acoustic  lenses  [57],  and  acoustic  rectifiers  [58], 

The  nonlinearity  of  the  interaction  law  results  from  the  Hertzian  contact  between  particles 
with  elliptical  contact  area  [6,  59,  60].  The  Hertzian  contact  relates  the  contact  force  Fy+i 
between  two  particles  (i  and  i+1)  to  the  relative  displacement  Ay+j  of  their  particle 
centers,  as  shown  in  the  following  equation 


Figure  1.  Granular  crystals  in  one-,  two-  and  three-dimensions  composed  by  metallic 
particles  confined  by  supporting  walls  or  confined  in  a  matrix. 


Values  inside  the  bracket  [s]+  only  take  positive  values,  which  denotes  the  tensionless 
characteristic  of  the  system  (i.e.,  there  is  no  force  between  the  particles  when  they  are 
separated).  For  Ay+^O  the  particles  are  just  touching,  Ali+i>0  the  particles  are  in 
compression,  and  Aii+i<0  the  particles  are  separated.  This  tensionless  characteristic  is  one 
part  of  the  nonlinearity  of  the  Hertzian  contact. 

For  two  spheres  (or  a  sphere  and  a  cylinder): 


,  (1) 

where  Et ,  vt ,  Rt  are  the  elastic  modulus,  the  Poisson's  ratio,  and  the  radius  of  the  i-th 

particle,  respectively.  The  ni  i+1  =3/2  comes  from  the  geometry  of  the  contact  between 

two  linearly  elastic  particles  with  elliptical  contact  area,  as  can  be  seen  in  [59].  In 
addition  to  assuming  the  contact  area  is  elliptical,  and  that  both  particles  remain  linearly 
elastic,  the  derivation  of  Hertzian  contact  assumes  [59]  (i)  the  contact  area  is  small 
compared  to  the  dimensions  of  the  particle,  (ii)  the  contact  surface  is  frictionless  with 
only  normal  forces  between  them,  (iii)  the  motion  between  the  particles  is  slow  enough 
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that  the  material  responds  quasi-statically.  Variation  of  the  contact  geometry  can  results 
in  variations  of  the  interaction  law  stiffness  and/or  nonlinearity,  and  ultimately  in 
variation  of  the  acoustic  properties  of  the  crystals.  Recent  works  studied  theoretically, 
numerically  and  experimentally  the  dynamic  response  of  chains  of  particles  composed  of 
grain  with  different  geometries  [61-63]  (see  section  6.3.4). 

The  rest  of  this  chapter  describes  past  and  recent  work  in  one-dimensional  (ID)  and  two- 
dimensional  (2D)  granular  crystals,  categorized  according  to  periodicity  and  dynamical 
regime. 


6.2  One-dimensional  Granular  Crystals. 

The  dynamic  properties  of  one-dimensional  (ID)  granular  crystals  have  been  extensively 
studied,  using  analytical,  numerical,  and  experimental  methods.  In  the  following  sections 
we  describe  some  of  the  most  interesting  physical  phenomena  supported  by  these 
nonlinear  systems. 

If  the  stiffness  of  the  contact  between  two  adjacent  particles  is  very  low  compared  to  the 
bulk  stiffness  of  the  particles  composing  the  crystal,  ID  granular  crystals  can  be  modeled 
as  a  system  of  nonlinear  springs  and  point  masses  (FPU-like  nonlinear  lattices).  Another 
perspective  from  which  to  approach  this  same  idea  is  that  the  characteristic  (resonant) 
frequencies  of  the  particles  themselves  are  very  high  compared  to  the  frequencies  of  the 
granular  crystal  system  involving  the  rigid  body-like  motion  of  the  particles  in  the 
system.  Neglecting  any  dissipation,  a  statically  compressed  ID  array  of  elastic  granules  is 
described  by  the  following  system  of  coupled  nonlinear  differential  equations: 


mMi  =  4K  +  ui-i  ~  R]+  -4+i Km  +  ui  ~  um~\T  •  (2) 


For  spherical  particles  we  recall  that  nj  =  —  and  Aj  is  defined  as  in  equation  (1).  Here,  the 
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static  overlap  S0i  =  —  ,  and  F0  is  the  homogeneous  static  compression  force.  mj  is 

’  U J 

the  mass  of  the  zth  particle  and  u,  is  the  dynamic  displacement  of  the  i  th  particle  from  its 


equilibrium  position  in  the  initially  statically  compressed  chain.  The  bracket  [s]+  of  Eq. 
(*)  takes  the  value  5  if  s  >  0  and  the  value  0  if  s  <  0,  which  signifies  that  adjacent  beads 
are  not  in  contact.  Within  this  framework,  the  dynamic  of  the  system  can  be  tuned  to 
encompass  linear,  weakly  nonlinear,  and  strongly  nonlinear  regimes  of  dynamic 
behavior,  as  will  be  introduced  for  the  mono-atomic  case  in  the  following  section. 


6.3  One-dimensional  Monoatomic  Granular  Crystals. 


We  start  by  studying  the  nonlinear  dynamic  behavior  of  a  statically  compressed  ID 
monoatomic  granular  crystal  (all  particles  are  the  same).  We  consider  the  case  of  a 
granular  crystal  composed  of  identical  elastic  spherical  granules,  as  shown  in  figure  2, 
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namely  Rt  -  R ,  nr  =  m  =  —  nR  p{) ,  and  A,  of  equation  (2)  is  reduced  to 


A  =  A  = 


EsjlR 

sp) 


,  where  m  is  the  mass  of  the  sphere,  E  and  po  are  the  Young’s  modules 


and  density  of  particle  material,  R  is  the  particle  radius  and  v  is  the  Poisson’s  ratio. 
Moreover,  we  assume  that  the  chain  is  subjected  to  constant  compression  forces  F0 
applied  to  both  ends,  resulting  in  an  initial  displacement  S0  between  neighboring  particle 


centers,  S0i  = 


r  \2/3 

F  \ 

—  =  S0  .  The  particle  equations  of  motion,  shown  in  equation  (2)  thus 

j 


reduce  to: 


mu,  =  A[S0  +  - u;]f  -  A[S0  +  ul  -  w,.+1]f  ,  (3) 

where  ui  is  the  displacement  of  the  i  th  bead  from  its  equilibrium  position  in  the  initially- 
compressed  chain,  as  shown  in  figure  2,  and  i  e  {2,  •  •  • , N  - 1} . 


F0  Uj-|  Ui  Uj  +■  i  F0 


Figure  2.  One-dimensional  monoatomic  crystal  compressed  by  a  static  force  F(j .  The 
crosses  represent  the  initial  positions  of  the  particles  centers  in  a  statically  compressed 
chain  while  the  black  circles  denote  the  current  positions.  The  right  end  of  the  chain  is 
undisturbed.  The  direction  of  the  impulse  propagation  is  shown  by  arrow  [13]. 


6.3.1  Near-linear  regime. 


To  approximate  the  fully  nonlinear  equations  of  motion  shown  in  equation  (3)  we  can 
take  a  power  series  expansion  of  the  forces.  For  dynamical  displacements  with  amplitude 

U  —  ll- 

much  less  than  the  static  overlap,  i.e.  — — - —  <sc  1 ,  one  can  keep  only  the  hannonic 

term  of  the  expansion.  In  this  case,  the  granular  crystal  can  be  considered  as  a  linear 

3 

lattice  with  spring  constants,  =  —A5q2  ,  where  the  equations  of  motion  are  reduced  to: 

miui=K2.(ui_l-Ui)-K2M(Ui-Ui+l) .  (4) 

The  spectral  band  of  the  ensuing  linear  chain  (see  chapter  2  for  more  details)  has  an  upper 
cutoff  frequency  of  com  =  yj4 K2  /  m  .  As  a  consequence  of  the  nonlinear  relation  F0 «  c)3 2 
for  the  case  of  the  spherical  granules,  the  cutoff  frequency  (as  well  as  the  sound  velocity 
of  the  ID  monoatomic  granular  crystal)  scales  as  /y  6 .  These  results  have  been  confirmed 
quite  accurate  experimentally  [46,  64]. 


6.3.2  Weakly  nonlinear  regime. 

\  ll  —  'll  \ 

If  the  dynamic  displacements  have  small  amplitudes  1  relative  to  those  due  to 

*o, 

static  load,  we  can  take  the  power  series  expansion  of  the  forces  (up  to  quartic 
displacement  terms)  to  yield  the  K2-  K3-  K4  model: 

mil,  =  K2(u :m  -  2m,.  +  u,_, )  +  K3  ((n,+1  -  n,.)2  -  (km  -  uif)  +  K4  ((ni+1  -  m,.)3  +  (t/M  -  m,)3)  ,(5) 
where  K2  =  ^Ad"\  K,  =  - K 4  =  ^AS~la . 

This  model  is  an  example  of  the  celebrated  FPU  model.  Many  theoretical  studies  have 
focused  on  the  dynamical  properties  of  this  type  of  nonlinear  lattice,  revealing  the 
existence  of  coherent  nonlinear  structures  such  as  nonlinear  periodic  waves,  solitary 
waves  [65],  and  discrete  breathers  [66], 

Seeking  traveling  waves  with  a  characteristic  spatial  size  L  that  is  much  larger  than  the 
inter-particle  distance  a  =  2R  -  §o,  one  can  apply  the  so  called  long-wavelength  or 
continuum  approximation.  Using  the  replacement: 


Ui  =  um  ~u±  aux  +  +  T7a\x .  (6) 
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equation  (5)  is  transformed  into  the  nonlinear  Boussinesq  equation  and  into  the 
Korteweg-de  Vries  equation  (see  for  example  [33]).  Nesterenko  applied  this  method 
(taking  into  account  only  the  K 3  term)  to  a  strongly  compressed  granular  chain,  and 
derived  the  following  KdV  equation: 

£ + co4v  +  r4m + -%-&x  =  o,  4  =  ~u, 

2c0 
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The  solutions  of  equation  (7)  are  well  known,  and  include  nonlinear  periodic  waves  and 
solitary  waves. 

On  the  other  hand,  by  investigating  how  quasi-monochromatic  plane  waves,  or  narrow- 
band  packets  evolve  by  nonlinear  effects,  one  can  derive  another  well-known  nonlinear 
wave  equation  -  the  Nonlinear  Schrodinger  (NLS)  equation.  This  equation  predicts  many 
nonlinear  phenomena,  including:  second  hannonic  generation,  modulation  instability,  and 
the  existence  of  bright  and  dark  solitons  [34].  The  derivation  of  the  NLS  from  equation 
(5)  is  possible  using  the  method  of  multiple  scales  combined  with  a  quasi-discreteness 
approximation  (see  [67]  for  a  generic  FPU  lattice  of  the  form  of  equation  (5)  and  [68]  for 
a  recent  application  of  this  method  to  a  monoatomic  strongly  compressed  granular 
crystal). 

Another  generic  feature  of  nonlinear  lattices  is  the  discrete  breather  (DB).  DBs  have  been 
studied  extensively  in  monoatomic  FPU  chains  [66].  One  of  the  mechanisms  for  the 
generation  of  such  nonlinear  localized  modes  is  the  modulational  instability  (MI)  of  the 
plane  wave  at  the  band  edge.  A  detailed  analysis  of  this  instability  (bifurcation)  shows 
that  the  MI  of  the  upper  cutoff  mode  manifests  itself  when  3 K2K4  -  4 K\  >  0  (see  Section 

4.3  of  Ref.  [31]  and  references  therein).  In  the  monoatomic  granular  crystal  setting,  one 
can  show  that  this  inequality  does  not  hold.  This  is  an  indication  that  small-amplitude 
DBs  bifurcating  from  the  upper  band  mode  do  not  exist  in  monoatomic  granular  crystals. 
However,  the  existence  of  dark  discrete  breathers  or  large-amplitude  DBs  remains  an 
interesting  open  question.  Finally,  one  more  interesting  weakly  nonlinear  effect,  the  self¬ 
demodulation,  was  studied  by  Toumat  and  collaborators  in  compressed  ID  granular 
crystals  [69]. 


6.3.3  Highly  nonlinear  regime:  long-wavelength  approach. 


A  very  interesting,  non-classical  wave  behavior  appears  if  the  granular  material  is  weakly 
compressed  and  the  neighboring  particle  displacements  are  larger  than  the  initial  relative 
displacement  do,  resulting  from  the  static  compression.  We  consider  this  regime,  the 
highly  nonlinear  regime.  Most  of  the  studies  in  ID  monoatomic  granular  crystals  have 
been  devoted  to  this  dynamical  regime.  In  this  section,  we  present  the  basic  steps  of  the 
long-wavelength  method  that  Nesterenko  applied.  A  review  of  alternate  analytical 
approaches  and  experimental  observations  will  be  presented  in  the  following  sections. 


Figure  3.  Weakly  compressed  chain  of  particles.  The  crosses  represent  the  initial 
positions  of  the  particles  in  the  statically  compressed  chain,  the  black  circles  correspond 
to  the  current  positions  of  spheres,  and  the  open  circles  the  initial  positions  of  the  spheres 
in  the  uncompressed  chain.  The  right  end  of  the  chain  is  undisturbed  and  the  direction  of 
the  impulse  propagation  is  shown  by  arrow  [13]. 


Including  So  in  displacement  uj  calculated  from  the  particle  positions  in  the 
uncompressed  system  (see  [13]  for  more  details),  equation  (3)  becomes: 

mu,  =  A[u,  ,  -  u,f'2  -  A[u,  -  uM ]f  (8) 

In  the  long-wavelength  approximation,  the  displacements  ui+j  can  be  expanded  in  a 
power  series  according  to  a  small  parameter  s  =  a  /L  up  to  the  fourth  order  (see  equation 
(6)).  By  substituting  (6)  into  equation  (8),  and  conducting  some  additional  calculation,  a 
new  wave  equation  is  obtained: 


Despite  the  complex  nature  of  the  presented  strongly  nonlinear  wave  equation,  the 
stationary  solutions  of  equation  (9),  such  as  nonlinear  periodic  and  solitary  waves,  can  be 


found  in  the  fonn  u(x  -  Vt)  [13].  The  wavefonn  of  a  periodic  wave  with  speed  Vp  is  given 
by  the  following  expression: 
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The  dependence  of  the  speed  of  the  periodic  wave,  Vp,  on  minimal  and  maximum  strains 
£  =  -ux  Cm.n,  Cmax)  is  presented  in  [13]. 


The  solitary  shape  (for  the  case  when  the  initial  prestrain  is  approaching  0)  is  one 
hump  of  the  periodic  solution  of  equation  (10),  with  finite  wave  length  equal  only  five 
particle  diameters.  This  solitary  wave  is  a  supersonic  one  as  well  as  a  KdV  soliton.  A 
unique  feature  of  this  solitary  wave  is  the  independence  of  its  width  on  amplitude. 
Accordingly,  this  property  is  quite  different  from  the  property  of  weakly  nonlinear  KdV 
solitary  wave.  Here,  the  speed  of  the  solitary  wave  Vs  has  a  nonlinear  dependence  on 
maximum  strain  2rn  (and  particle  velocity-nm): 
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This  result  shows  that  the  speed  of  the  strongly  nonlinear  solitary  wave  Vs  does  not 
depend  on  particle  size  in  the  granular  material.  At  the  same  time  it  does  depend  on  the 
elastic  properties  of  the  particles  ( E  and  v)  and  their  density.  The  presented  theoretical 
results  allow  us  to  design  strongly  nonlinear  granular  materials  with  exceptionally  low 
velocity  of  signal  propagation.  Simple  estimation  based  on  equation  (11)  shows  that  it  is 
possible  to  create  materials  with  nonlinear  impulse  speed  in  the  interval  10  -  100  m/s 
corresponding  to  audible  signal  speeds. 


6.3.4  Review  of  alternate  strongly-nonlinear  wave  theoretical  approaches. 

The  presented  solitary  wave  solution  {soliton  with  compact  support  known  also  as 
compacton  [70])  describes  very  well  the  solitary  wave  that  a  impulsive  excitation 
generates  in  a  weakly  compressed  or  uncompressed  chain.  This  was  verified  in 
simulations  and  experiments  by  different  authors  (see  references  below).  The  rigorous 
proof  of  the  existence  of  solitary  waves  in  a  monoatomic  granular  crystal  composed  of 
spherical  particles  was  done  by  MacKay  [71]  ,  who  applied  the  existence  theorem  for 
solitary  waves  on  lattices  by  Friesecke  and  Wattis  [72].  Ji  and  Hong  extended  the  proof 
given  by  MacKay  to  the  general  case  of  an  arbitrary  power-law  type  contact  force  [73]. 

An  analytical  solution  of  the  form  tanh(fB )  for  stationary  waves  in  discrete  chains,  where 
fD  is  represented  by  a  series,  is  presented  by  Sen  and  Manciu  [74].  Their  result  is  very 


close  to  a  soliton  obtained  by  the  long-wavelength  approximation.  Chatterje  studied  the 
asymptotic  description  of  the  tail  of  the  soliton  in  an  uncompressed  chain  and  he  revealed 
its  double  exponential  decay  [75].  He  also  presented  a  new  asymptotic  solution  for  the 
full  solitary  wave,  which  is  closer  to  the  results  of  numerical  simulations  than  the 
approximate  solution  given  by  Nesterenko.  A  quite  different  analytic  approach  for  the 
study  of  pulse  propagation  in  granular  crystals  was  developed  by  Lindenberg  and 
collaborators  [76,  77].  This  method  uses  the  binary  collision  approximation  to  reduce  the 
problem  of  propagation  to  collisions  involving  only  two  granules  at  a  time. 

English  and  Pego  [78]  studied  the  shape  of  the  solitary  wave  that  propagates  in  a  Id 
granular  chain  without  precompression  (8o=0).  Their  method  is  based  on  a  reformulation 
of  the  equations  of  motion  using  the  difference  coordinates  r  =  ujl  -  uj  such  that: 

mri=A[[rM]f-2[ri]f  +  [ri^]f\  (12) 

Seeking  for  traveling  wave  solutions,  r  =  r(x)  =  r(i  -  ct )  ,  one  obtains  the  following 
advanced  delay  equation: 

r"(x )  =  — ^-Tr3/2(x  - 1)  -  2 rV2(x)  +  rv\x  + 1)1 .  (13) 
me  L  J 

By  rewriting  this  equation  in  an  equivalent  integral  form  and  studying  its  asymptotic 
behavior  they  proved  that  the  solitary  wave  decays  super-exponentially.  Moreover,  they 
applied  an  iterative  method  for  the  computation  of  the  numerically  exact  shape  of  the 
solitary  waves. 

Later,  Ahnert  and  Pikovksy  [79]  applied  a  different  type  of  quasi-continuum 
approximation  by  expanding,  up  to  fourth  order,  the  difference  coordinate  r  instead  of 

the  displacement  uj .  Substituting  these  expansions  in  equation  (12),  they  obtained  a 

strongly  nonlinear  partial  differential  equation  (see  equation  (6)  in  [79])  that  supports  a 
solitary  wave  with  compact  fonn.  The  analytic  solution  has  the  same  form  with 
Nesterenko’s  solution,  but  with  slightly  different  amplitude  and  width  constants. 
Moreover,  they  presented  an  accurate  numerical  method  for  the  numerical  solution  of  the 
advanced  delay  equation  (13)  and  they  compared  the  numerically  obtained  solutions  with 
those  of  approximated  PDEs. 

Recently,  Starosvetsky  and  Vakakis  [80],  working  directly  on  the  nonlinear  lattice 
equations  with  no  precompression,  developed  semi-analytical  approaches  for  computing 
different  families  of  nonlinear  traveling  waves.  These  waves  involve  both  separation  and 
compression  between  adjacent  particles  and  therefore  they  cannot  be  resolved  using 
quasi-continuum  approximations.  In  addition,  they  showed  that  these  wave  families 
converge  to  the  solitary  wave  in  a  certain  asymptotic  limit.  They  also  solved  the  reduced 
advanced  delay  equation  numerically,  and  applied  the  method  of  Pade  approximations. 


6.3.5  Review  of  experiments  with  strongly  nonlinear  solitary  waves. 

The  quantitative  agreement  of  analytical  and  numerical  predictions  with  experiments  on 
strongly  nonlinear  granular  crystals  was  first  found  by  Lazaridi  and  Nesterenko,  [44]. 
They  observed  for  the  first  time  the  rapid  decomposition  of  the  initial  impulse  excitation 
into  multiple  solitary  waves  in  a  distance  comparable  to  the  solitary  wave  width.  Optical 
observations  of  strongly  nonlinear  waves  were  reported  by  Zhu,  Shukla  and  Sadd  [81]. 
They  reported  the  fonnation  of  nonlinear  waves  in  photoelastic  disks  excited  by  a  local 
explosive  loading.  Coste,  Falcon  and  Fauve,  [45],  Coste  and  Gilles,  [46]  conducted  a  very 
detailed  quantitative  study  of  the  speed  and  shape  of  solitary  waves  at  different  solitary 
wave  amplitude.  They  reported  a  negligible  decay  of  the  solitary  wave  in  chains 
composed  of  50  particles,  and  they  concluded  that  the  solitary  waves  shape  observed  in 
experiments  is  in  very  good  agreement  with  the  predictions  obtained  from  theoretical 
solutions,  equation  (10). 
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Figure  4.  Evolution  of  the  soliton  train  excited  in  experiments  by  striker  impact  (Ms=10/u, 
us=0.5  m/s)  with  the  propagation  distance:  (a)  N  =  5,  (b)  N  =  10,  (c)  N  =  20,  (d)  N  =  30, 
(e)  N  =  40,  (f  )  N  =  60.  Vertical  scale  corresponds  to  80  Newton,  horizontal  scale  to  50 
ps  (a),(b),(c),(d),(e)  and  100  ps  (f),  N  is  number  of  particles  [13]. 


It  should  be  emphasized  that  the  relatively  low  speed  of  the  waves  (500  m/s)  detected  by 
Coste  et  ah,  [45]  is  very  unusual  for  solid  materials  and  corresponds  to  a  disturbance  of 
relatively  large  amplitude.  Theoretically  the  minimum  propagation  speed  of  a  solitary 
wave  can  be  close  to  zero  if  the  amplitude  of  the  disturbance  is  approaching  zero, 
equation  (11).  Using  polymeric  and  composite  particles,  for  example,  one  can  design 
granular  crystals  with  a  solitary  wave  speed  corresponding  to  an  audible  signal  in  the 
interval  of  10  -  100  m/s,  an  order  of  magnitude  less  than  previously  observed  in 
experiments  by  Coste  et  ah,  1997  [82].  Nesterenko  et  al.  showed  experimentally  the 
presence  of  anomalous  reflections  when  highly  nonlinear  waves  interact  with  interfaces 
[52],  effectively  demonstrating  for  the  first  time  the  concept  of  an  acoustic  diode.  Daraio 
et  al.,  2006  [47]  described  in  detail  the  ability  to  tune  the  dynamic  response  of  granular 
crystals  by  controlling  the  static  precompression  and  the  dynamic  excitation  applied  to 
the  system.  Job  et  al.,  2005,  [83]  investigated  the  behavior  of  solitary  waves  interacting 


with  a  boundary,  showing  for  the  first  time  the  sensitivity  of  solitary  waves  to  the 
mechanical  properties  of  an  adjacent  medium.  A  thorough  experimental,  numerical  and 
theoretical  descriptions  of  the  formation,  propagation  [84],  and  collision  of  solitary  wave 
trains  were  published  a  few  years  later  [85], 

Recently,  several  experimental  works  described  the  wave  propagation  in  granular  crystals 
composed  of  particles  with  elliptical  and  cylindrical  geometry  [61,  62],  For  example, 
chains  composed  of  ellipsoidal  or  cylindrical  particles  were  shown  to  support  the 
fonnation  and  propagation  of  highly  nonlinear  solitary  waves  similar  to  the  solitary 
waves  observed  in  chains  of  spherical  particles.  However,  these  systems  showed  to  be 
also  highly  dependent  on  the  particles  geometry  and  on  the  orientation  angles  between 
particles  in  the  chain.  This  interesting  dependence  on  orientation  angle  between  beads 
provides  an  additional  free  parameter  to  design  acoustic  materials  with  unprecedented 
transmission  properties. 

Experimental  data  obtained  also  tudying  the  dynamic  response  of  a  chain  composed  of 
spherical  steel  particles  coated  with  a  soft  polymeric  material  showed  that  these  systems 
also  support  the  fonnation  and  propagation  of  highly  nonlinear  solitary  waves  [82]. 
However,  one  interesting  property  of  these  systems  is  that  the  contact  interaction  between 
thin-coated  spheres  does  not  follow  classical  Hertzian  interaction  between  two  solid 
spheres  [86].  The  dynamic  response  of  chains  of  coated  spheres  is  governed  by  a 
quadratic  power  law  dependence  between  the  contact  force,  F,  and  the  displacement,  6, 
instead  of  the  Hertzian,  non-integer  power  of  3/2.  This  new  nonlinear  contact  interaction 
dramatically  changes  the  dynamics  of  solitary  wave  propagation  compared  to  its 
counterpart  in  chains  of  solid  spheres.  Here,  the  spatial  width  of  the  wave  becomes 
shorter  (3.14  particles  size  instead  of  5)  the  wave  speed  (Vs)  is  relatively  slow,  and  its 
dependence  on  force  amplitude  (Fm)  is  also  different  (Vs  ~  Fm  1/4  instead  of  Vs  ~  Fm  1/6). 

Similarly  studies  of  chains  of  hollow  spherical  particles  presented  interesting  nonlinear 
acoustic  phenomena.  Highly  nonlinear  solitary  waves  were  observed  to  propagate 
through  the  system,  but  the  wave  properties  were  different  from  the  highly  solitary  waves 
in  the  chains  of  solid  spheres.  The  spatial  width  of  the  solitary  wave  in  chain  of  hollow 
sphere  was  approximately  8  particles  size  (larger  than  5  beads  size),  and  the  wave  speed 
was  proportional  to  force  amplitude  to  power  1/1 1  [63].  It  was  shown  that  such  behavior 
derived  from  the  unique  contact  interaction  between  thin  hollow  spheres,  which  for  the 
range  of  wave  amplitude  studied  could  be  approximated  by  a  power-law  type  relation 
(. F=kdn ).  In  this  case,  the  exponent  n  was  found  to  be  smaller  than  the  value  3/2  as  in  the 
classical  Hertzian  interaction  between  solid  spheres.  The  contact  stiffness  k  and  the 
exponent  n  were  also  found  to  be  dependent  on  the  thickness  of  the  hollow  sphere’s  shell. 
This  dependence  of  the  dynamic  behavior  of  granular  crystals  on  the  coating  and/or  shell 
thickness  of  spherical  particles  provides  yet  another  free  parameter  to  employ  in  tuning 
the  dynamics  of  nonlinear  acoustic  crystals. 


6.4  One-dimensional  Diatomic  Granular  Crystals. 


By  increasing  the  degree  of  periodicity,  from  a  homogenous  monoatomic  granular  crystal 
to  a  diatomic  granular  crystal  composed  of  alternating  particles,  we  gain  access  to 
additional  interesting  phenomena.  This  section  describes  some  of  those  phenomena 
characteristic  of  ID  diatomic  granular  crystals,  including  tunable  band  gaps,  discrete 
(DBs),  and  highly  nonlinear  solitary  waves  with  widths  up  to  10  particles. 

An  example  of  a  ID  diatomic  granular  crystal  is  illustrated  in  the  bottom  of  figure  1. 
Here  we  reduce  the  equation  of  motion  for  the  general  ID  granular  crystal,  shown  in 
equation  (2)  ,  to  the  ID  diatomic  crystal  model,  as  follows: 

mu  =  A[S0  +  -  w,]f  -  A[S0  +  ui  - uMr ,  (14) 

where  the  subscript  i  is  the  index  of  the  7th  particle,  the  particle  masses  are  m2M  =  in  and 
m2j  =  M  .  By  convention,  we  will  take  M  to  be  the  larger  of  the  two  masses  and  m  to  be 
the  smaller  of  the  two  masses.  Because  all  contacts  (aside  from  any  boundaries)  are  the 
same,  there  is  a  single  Hertzian  contact  coefficient  A  and  static  overlap  So  that  are  used  to 
represent  the  system,  which  have  been  defined  in  the  previous  sections.  Within  this 
framework,  as  before,  the  dynamic  of  the  system  can  be  tuned  to  encompass  linear, 
weakly  nonlinear,  and  strongly  nonlinear  regimes  of  dynamic  behavior.  Also  as  before, 
the  K2-K3-K4  model  can  be  applied  in  the  weakly  nonlinear  regime,  and  the  AT  linearized 
model  in  the  linear  regime. 


6.4.1  Near-linear  regime:  localized  surface  modes 

For  dynamical  displacements  with  amplitude  much  less  than  the  static  overlap 
(|  uM  -  ut  |«:  S0) ,  we  can  neglect  the  nonlinear  Kx  and  K4,  and  compute  the  linear 

dispersion  relation  of  the  system.  This  results  in  an  effectively  linear  diatomic  system  of 
springs  and  point  masses,  as  was  presented  in  chapter  2,  but  with  a  tunable  stiffness  Ki. 

Several  previous  studies  explored  the  existence  of  band  gaps  in  highly  compressed 
granular  crystals.  Initially,  studies  focused  on  ID,  with  a  two  particle  unit  cell,  arrays  of 
glued  [87],  welded  [88],  and  elastically  compressed  spherical  particles  [48,  54,  89].  These 
studies  demonstrated  tunable  vibration  spectra  with  two  bands  of  propagation  (called  the 
acoustic  and  optical  bands)  separated  by  a  band  gap  in  the  diatomic  case.  Boechler  et  al. 
[49]  later  extended  this  work  by  investigating  the  response  of  one-dimensional  diatomic 
granular  crystals  with  three-particle  unit  cells,  and  showing  their  tunability  based  on 
variations  of  the  particles  geometry  and  on  the  applied  static  load.  In  contrast  to  diatomic 
granular  crystal  with  two  particle  unit  cells,  the  three  particle  unit  cell  granular  crystal 
was  shown  to  contain  up  to  three  distinct  pass  bands  and  a  two  finite  band  gaps. 


In  addition  to  acoustic  and  optical  band  modes,  the  diatomic  semi-infinite  hannonic 
granular  crystal  also  supports  a  gap  mode,  provided  the  existence  of  a  light  particle  at  the 
surface  and  free  boundary  conditions.  This  mode  is  localized  at  the  surface  (i.e.,  at  the 
first  particle)  and  its  displacements  have  the  following  fonn  [90]: 
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with  particle  number  /'  >  0 ,  frequency  cos  =  ^ K2(\  /  in  +  1  /  M)  is  in  the  gap  of  the  linear 


spectrum,  and  B  is  an  arbitrary  constant.  This  particular  mode  with  frequency  in  the  band 
gap,  localized  around  the  surface,  proves  to  have  a  nonlinear  counterpart  and  to  be  very 
closely  related  to  the  DB  in  the  strongly  discrete  regime  as  will  be  described  in  the 
following  section. 


6.4.2  Weakly  nonlinear  regime:  discrete  breathers 

By  increasing  the  relative  amplitude  of  the  dynamic  to  static  displacements 
(\uM-ui\<S0) ,  and  thus  increasing  the  nonlinearity  of  the  response  into  the  weakly 

nonlinear  regime,  a  type  of  intrinsic  localized  mode  called  a  discrete  breather  (DB)  can 
be  supported  by  the  system.  DBs  have  been  widely  studied  in  the  realm  of  nonlinear 
lattices  as  previously  described  [31].  They  are  nonlinear  modes  that  have  frequency 
within  the  gap  of  the  linear  spectrum  and  are  localized  in  space.  As  such,  discrete 
breathers  have  practical  importance  as  a  mechanism  to  localize  vibrational  energy  in 
frequency  and  space  without  the  introduction  of  any  extrinsic  disorder. 

DBs  were  rigorously  proven  to  exist  in  diatomic  FPU-type  lattices,  with  alternating  heavy 
and  light  masses,  by  Mackay  in  1997  [91],  Furthermore,  several  studies  also  investigated 
the  specific  case  of  DBs  located  in  the  gap  between  the  acoustic  and  optical  bands  of 
anharmonic  diatomic  lattices  [92-94]. 

A  recent  study  by  Theocharis  et.  al.  systematically  studied  the  existence  and  stability  of 
DBs  in  diatomic  granular  crystals  [51],  Studies  in  other  diatomic  anhannonic  lattices 
have  shown  the  existence  of  up  to  two  types  of  DBs.  The  study  by  Theocharis  et.  al. 
demonstrated  that  both  can  arise  in  granular  chains.  They  examined  both  of  these  two 
families  of  discrete  gap  breathers,  and  studied  their  existence,  stability,  and  structure 
throughout  the  gap  of  the  linear  spectrum.  The  first  family  was  an  unstable  DB  that  is 
centered  on  a  heavy  particle  and  characterized  by  a  symmetric  spatial  energy  profile,  and 


the  second  family  is  a  potentially  stable  DB  that  is  centered  on  a  light  particle,  and  is 
characterized  by  an  asymmetric  spatial  energy  profile. 

Although  the  FPU  and  granular  crystal  lattices  are  analogous  in  many  respects,  there 
exists  an  important  difference  because  of  the  additional  nonlinearity  caused  by  the 
tensionless  characteristic  of  the  granular  crystal  lattice.  Accordingly,  Theocharis,  et.  al. 
contrasted  discrete  breathers  in  anharmonic  FPU-type  diatomic  chains  with  those  in 
diatomic  granular  crystals,  and  found  that  for  the  case  when  the  DB  was  very  narrow 
(highly  discrete)  that  the  asymmetric  nature  of  the  latter  interaction  potential  led  to  a  form 
of  hybrid  bulk-surface  localized  solutions  (see  figure  5).  Figure  5  shows  the  two  families 
of  DB  solutions  at  times  t  =  7  and  t  =  7/2  (where  7  is  the  periodic  of  the  DB)  with  the 
profile  of  a  linear  surface  mode.  This  similarity  between  the  shapes  of  the  two  modes 
suggests  that  the  temporary  creation  of  a  new  interior  surface,  allowed  by  tensionless 
characteristic  of  the  system,  has  contributed  to  a  modified  type  of  intrinsic  localized 
mode. 
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Figure  5:  Top  panels:  Spatial  profile  of  a  DB  in  the  heavy  mass  centered  symmetric 
family  a  times  (a)  t  =  0  and  (b)  t  =  7/2.  Bottom  panels:  As  with  the  top  panels,  but  for  the 
light  mass  centered  asymmetric  family  of  DB  solutions.  The  dashed  curves  correspond  to 
the  spatial  profile  of  the  surface  mode  obtained  using  equations  (15).  In  each  panel,  we 
include  a  visualization  of  particle  positions,  and  the  corresponding  spatial  gap  openings, 
for  the  corresponding  time  and  DB  solution. 


The  existence  of  DBs  in  diatomic  granular  crystals  was  experimentally  proven  in  a  recent 
study  by  Boechler,  et.  A1  [50].  In  this  study,  the  authors  utilized  the  modulational 
instability  (MI)  of  the  lower  optical  mode  to  generate  DBs  in  at  80  particle  diatomic 
granular  crystal.  In  the  weakly  nonlinear  regime,  granular  crystals  can  be  showed  to  be 
subject  to  MI  when  K3  /K.2K4<3/4.  To  excite  the  MI,  they  drove  the  granular  crystal  from 
one  boundary  at  the  lower  optical  mode  frequency,  at  high  amplitude.  Upon  reaching  a 
critical  amplitude  for  the  MI  to  occur,  the  harmonic  lattice  vibration  decayed  into  a 
localized  DB. 

In  figure  6,  as  per  Boechler,  et.  al.  [50],  we  show  an  experimental  observation  of  a  DB, 
generated  in  an  80  particle  diatomic  granular  crystal.  This  example  shows  how  the 
interplay  of  nonlinearity  and  discreteness/periodicity  leads  to  the  localization  of 
vibrational  energy  to  a  narrow  spatial  regime  (around  the  14th  particle  from  the  boundary) 
at  a  specific  frequency  within  the  gap  of  the  linear  spectrum  (/y,  =  8.31  kHz).  Here  we  can 
see  in  panels  (a)  and  (b),  far  from  the  center  of  the  DB,  a  periodic  response  at  the  driving 
frequency  (fd=  8.9  kHz).  Alternatively,  in  panels  (c)  and  (d),  near  the  center  of  the  DB, 
we  see  a  quasiperiodic  response  characterized  by  the  driving  frequency  and  the  secondary 
DB  frequency.  This  spatial  localization  is  further  clarified  in  the  spatial  profile  of  the 
energy  distribution  shown  in  panel  (e). 
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Figure  6:  Experimental  observation  of  a  DB,  in  an  80  particle  granular  crystal,  at  fb  =  8.31 
kHz.  (a),(b)  Force  at  particle  2  and  14,  respectively.  (c),(d)  Power  spectral  density  (PSD) 
for  the  highlighted  time  regions  in  (a),(c)  of  the  same  color.  Square  [circular]  markers 
denote  the  DB  [driving]  frequency  and  PSD  amplitude,  (e)  The  ratio  of  the  PSD 
amplitude  at  the  discrete  breather  frequency  divided  by  the  PSD  amplitude  of  the  driving 
frequency  as  a  function  of  sensor  location.  The  vertical  dashed  line  in  (b)  and  (d)  denotes 
the  lower  cutoff  frequency  of  the  optical  band,  and  the  vertical  dashed  lines  in  (a)  and  (c) 
denote  the  time  region  for  the  PSD  calculation. 


6.4.3  Highly  nonlinear  regime:  strongly  nonlinear  solitary  waves 


In  this  section,  we  further  explore  the  effects  of  increased  periodicity  by  studying  the 
propagation  of  highly  nonlinear  solitary  waves  in  ID  diatomic  granular  crystals  with  no 
static  load.  Solitary  waves  in  such  systems  were  first  studied  and  described  by  V. 
Nesterenko  in  2001  [13].  He  found  that  by  assuming  the  mass  of  one  particle  type  is 
much  larger  than  the  mass  of  the  other  ( mi/m.2  »  1)  and  by  applying  the  long  wavelength 
approximation,  that  the  resulting  wave  equation  supports  a  solitary  wave  solution  with  a 
characteristic  spatial  width  of  ~10  particles.  This  demonstrates  how  an  increase  in 
periodicity  (or  redistribution  of  the  monoatomic  particle  masses  to  two  neighboring 
particles)  can  result  in  wider  solitary  wave. 

Later,  Porter  et  al.  [95]  applied  the  long-wavelength  approximation  to  diatomic  granular 
crystals  with  arbitrary  mass  ratios  by  postulating  a  “consistency  condition”  between  the 
displacements  of  the  two  particles  in  the  unit  cell.  They  showed  that  the  diatomic  chain 
supports  a  finite-width  soliton-like  solution,  and  they  obtained  an  analytical  expression 
for  the  width  of  the  solution  as  a  function  on  the  mass  ratio.  This  expression  generalizes 
the  previously  known  limiting  cases,  namely,  mj/m 2  =  1  (monoatomic)  and  m i/m?  »  1 
(diatomic  with  ~10  particle  length  solitary  wave  width).  In  the  same  study,  Porter  et.  al. 
compared  these  analytical  predictions  with  simulations  and  experiments  and  found  good 
agreement. 

Recently,  Vakakis  et.  al.  [96]  presented  an  extensive  numerical  and  theoretical  study  of 
solitary  waves  in  diatomic  chains.  They  showed  that  in  a  diatomic  granular  crystal, 
scattering  at  the  interfaces  of  the  dissimilar  light  and  heavy  beads  will  typically  cause  a 
slow  disintegration  of  the  traveling  wave  and  the  formation  of  small  amplitude  oscillating 
tails.  However,  they  also  found  that  for  specific  discrete  values  of  the  mass  ratio  between 
heavy  and  light  particles,  the  system  supports  solitary  waves  which  travel  without 
distortion.  These  discrete  values  of  the  mass  ratio  correspond  to  the  case  where  the  light 
beads  always  stay  in  contact  with  adjacent  heavy  beads.  For  this  case,  the  entire  energy  of 
the  main  pulse  is  conserved  and  transferred  without  loss  to  the  next  heavy  bead.  These 
solutions  can  be  considered  analogous  to  the  propagation  of  solitary  waves  in 
monoatomic  granular  crystals,  in  that  their  velocity  profiles  decay  to  zero.  Finally,  they 
also  observed  that  the  diatomic  family  of  solitary  waves  propagates  faster  than  the 
corresponding  solitary  waves  in  monoatomic  systems. 


6.5  One-dimensional  Monoatomic  Granular  Crystals  With  Defects 

By  introducing  one  or  more  defects  into  an  otherwise  perfectly  periodic  mono-atomic 
granular  crystal,  we  introduce  disorder  into  the  system.  The  effects  of  introducing 
disorder  into  the  system  are  analogous  to  the  effects  seen  in  the  previous  section,  as  the 


degree  of  periodicity  increased  from  monoatomic  to  diatomic  granular  crystals.  The 
presence  of  disorder,  and  its  interplay  with  the  nonlinearity  of  the  system,  causes  the 
existence  of  new  interesting  and  useful  phenomena  throughout  the  granular  crystal’s 
range  of  dynamic  regimes.  However,  in  contrast  to  the  case  of  increasing  periodicity,  by 
introducing  disorder  we  add  new  ways  to  break  the  spatial-symmetry  of  our  system.  In 
combination  with  the  ability  of  nonlinear  systems  to  break  the  time-reversal  symmetry  of 
the  dynamic  response,  the  introduction  of  spatially-asymmetric  disorder  can  be 
particularly  useful.  In  the  following  section  we  will  describe  several  recent  studies 
relating  to  defects  in  monoatomic  granular  crystals,  including:  tunability  of  defect  modes 
in  the  linear  regime  [97],  localized  nonlinear  defect  modes  and  spontaneous  symmetry 
breaking  in  the  weakly-nonlinear  regime  [98],  the  interplay  of  solitary  waves  with  defects 
in  the  highly  nonlinear  regime  [99],  and  tunable  bifurcation-based  acoustic  rectification 
in  a  driven  granular  crystal  [58]. 


6.5.1  Near-linear  regime:  tunable  defect  modes 

A  strongly-compressed  (with  respect  to  the  dynamic  displacements)  homogenous 
granular  crystal  with  light-mass  defects  will  contain  exponentially  localized  modes  with 
frequencies  above  the  acoustic  band  of  the  granular  crystal,  localized  around  the  defect 
sites.  The  frequency  of  these  localized  defect  modes  is  tunable  with  changes  in  static 
load,  similar  to  the  tunability  of  the  linear  dispersion  relation  of  a  periodic  granular 
crystal. 

The  existence  and  tunability  of  defect  modes  localized  around  one  and  two  light-mass 
defects  in  a  strongly-compressed  ID  otherwise  homogenous  granular  crystal  was 
investigated  first  numerically  and  analytically  by  Theocharis,  et.  al.  [51],  and  then 
experimentally  by  Man,  et.  al.  [97].  In  the  work  by  Man,  et.  al.,  they  placed  one  and  two 
light-mass  defects  near  the  edge  of  a  20  stainless-steel  particle  granular  crystal,  applied 
white-noise  excitation  from  the  edge  of  the  crystal,  and  measured  the  frequency  of  the 
defect  modes  localized  in  the  vicinity  of  the  defects  as  a  function  of  defect  size  and 
relative  defect  position.  The  observed  defect  mode  frequencies  were  compared  with 
eigen-analysis  of  the  linearized  20  particle  granular  crystal  (as  described  in  Theocharis, 
et.  al.),  and  analytical  expressions  based  on  few-site  considerations.  For  a  sufficiently 
small  single  light  mass  defect  in  an  otherwise  homogenous  granular  crystal,  the 
frequency  of  the  defect  mode  can  be  approximated  as  [97]: 
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which  is  obtained  by  solving  the  eigenvalue  problem  of  the  three-particle  system  in  the 
vicinity  of  the  defect  (large  particle-defect  particle-large  particle).  Here  M  is  the  mass  of 
the  homogenous  particles,  m  is  the  mass  of  the  defect  particle,  f?bead  is  the  frequency  of 
the  localized  defect  mode,  and  Krr=3/2Arr  Fo  '  is  the  linearized  stiffness  between  two 
large  particles  and  KRr=3/2ARr  Fo  is  the  linear  stiffness  of  the  contact  between  a 
defect-particle  and  a  large  particle.  ARr  and  Arr  are  the  Hertz  contact  coefficients 
between  the  respective  particles.  From  this  expression,  it  is  clear  how  the  defect  modes 
are  tunable  with  static  load,  geometry,  and  material  properties. 

In  both  studies  [97,  98],  it  was  found  that  when  two  defects  were  placed  sufficiently  far 
from  to  each  other  (outside  the  localization  length  of  each  individual  defect  mode),  that 
the  granular  crystal  presented  two  isolated  linear  defect  modes  with  frequencies  of  a 
single-defect  mode.  The  further  the  distance  between  the  defects,  the  closer  the  modes  are 
to  isolated  ones  with  near-identical  frequencies.  However,  when  the  defects  are  brought 
sufficiently  close  together  (within  the  localization  length  of  a  single-defect  mode)  each 
defect  was  found  to  affect  the  other.  This  caused  the  fonnation  of  a  symmetric  and  anti¬ 
symmetric  pair  of  defect  modes,  with  two  new  separate  frequencies,  involving  both 
defects. 


6.5.2  Weakly  nonlinear  regime:  nonlinear  localized  modes  and  symmetry  breaking 

If  we  increase  the  amplitude  of  the  dynamic  displacements,  relative  to  the  static  overlap, 
and  thus  increase  the  nonlinearity  of  the  dynamic  response,  the  nonlinear  localized  defect 
modes  depart  from  their  linear  counterparts  and  new  phenomena  are  introduced.  In 
addition  to  exploring  the  near-linear  behavior  of  one-dimensional,  strongly-compressed 
granular  crystals  with  one  or  two  light-mass  defects,  Theocharis,  et.  al.  investigated  the 
behavior  of  defects  in  the  weakly  nonlinear  regime  [98].  As  previously  described,  by 
analyzing  the  problem's  linear  limit,  they  identified  the  system  eigenfrequencies  and  the 
linear  defect  modes.  Using  continuation  techniques,  they  found  localized  nonlinear  defect 
mode  solutions  that  bifurcate  from  their  linear  counterparts  and  studied  their  linear 
stability  in  detail  by  computing  the  Floquet  multipliers  of  the  nonlinear  periodic 
solutions. 

For  the  case  of  a  single  light-mass  defect,  it  was  found  that  the  inherent  nonlinearity  of 
the  system  leads  to  long-lived  localized  breathing  oscillations,  which  fonn  robust 
nonlinear  localized  modes.  Their  frequency  depends  not  only  on  the  static  load,  the 
geometry  and  the  material  properties  of  the  granular  crystal  and  defect  particle,  but  also 
on  the  amplitude  of  the  oscillations.  Because  of  the  type  of  the  nonlinearity  in  the  system, 
the  defect  mode’s  frequency  decreases  with  increasing  dynamic  amplitude  (and 
nonlinearity).  Accordingly,  we  have  shown  two  ways  to  use  nonlinearity  to  tune  the 


frequency  of  a  localized  mode:  change  the  static  load  or  change  the  relative  amplitude  of 
the  dynamic  displacements  and  the  static  overlap. 

For  the  case  of  two  defects,  nonlinearity  can  create  further  interesting  phenomenology 
when  the  defects  are  sufficiently  close.  A  particularly  intriguing  example  is  the  case  of 
next-nearest  neighbor  defects,  where  the  two  defects  are  separated  by  one  large  particle. 
This  resembles  the  situation  of  a  “double  well”  potential,  which  has  been  studied 
systematically  in  various  settings,  including  nonlinear  optics  [100]  and  atomic  physics 
[101,  102].  In  these  settings,  it  has  been  predicted  analytically  (via  a  two-mode 
reduction),  manifested  numerically,  and  observed  experimentally  that  beyond  a  certain 
nonlinearity  threshold,  a  pitchfork  bifurcation  arises  that  causes  the  spontaneous 
symmetry  breaking  of  the  relevant  configurations,  and  results  in  the  predominant  support 
one  out  of  the  two  nonlinear  modes.  The  investigations  by  Theocharis,  et.  ah,  in  granular 
crystals,  indicate  that  this  phenomenology  is  generic.  Figure  7,  shows  the  bifurcation  of 
the  antisymmetric  linear  defect  mode,  as  a  function  of  the  defect  mode  frequency  and 
relative  force  felt  between  the  defect  sites,  for  the  a  next-nearest  neighbor  configuration. 
As  the  antisymmetric  defect  mode  (figure  7,  inset,  Al)  becomes  progressively  more 
nonlinear  (and  decreases  in  frequency),  at  a  critical  point  the  mode  becomes  unstable  via 
a  pitchfork-like  bifurcation.  The  bifurcation  signals  the  emergence  of  two  asymmetric 
modes  (figure  7,  insets  A2  and  A3),  which  are  mirror  images  of  each  other,  and 
predominantly  centered  on  one  of  the  two  defect  sites. 


Figure  7:  Pitchfork  bifurcation  illustrated  by  force  differential  between  two  next-nearest 
neighbor  defects,  as  a  function  of  the  mode  frequency.  This  shows  the  transition  from  a 
single  antisymmetric  mode  to  two  (mirror-symmetric  between  them)  asymmetric  modes 
after  the  onset  of  the  symmetry  breaking  bifurcation.  Insets:  spatial  profiles  and  locations 
of  Floquet  multipliers  A.  in  the  complex  plane  of  solutions  for  different  frequencies. 


The  case  of  the  bifurcation  of  the  antisymmetric  two-defect  mode  is  a  good  example  of 
how  through  the  addition  of  nonlinearity,  we  can  create  sharp-transitions  between  two 
acutely  different  states,  break  the  spatial  symmetry  of  the  dynamic  response,  and  access 
new  mechanisms  to  control  the  distribution  and  frequency  of  vibrational  energy. 


6.5.3  Highly  nonlinear  regime:  transient  localized  modes 

By  increasing  the  nonlinearity  of  the  dynamic  response  further,  we  can  investigate  how 
the  interaction  of  traveling  waves  with  defects  in  a  nonlinear  system.  The  interaction  of 
highly  nonlinear  solitary  waves  with  a  mass  defect  placed  in  a  ID,  unloaded  granular 
crystal  has  been  investigated  analytically  and  computationally  in  Hascoet,  2000  and 
numerically  and  experimentally  in  Job,  2009.  Two  different  physical  pictures  emerge 
whether  one  considers  a  light  or  a  heavy  impurity  mass.  The  scatter  of  the  solitary  wave 
with  a  light  impurity  yields  transient  oscillations  of  the  defect  which  leads  to  the  emission 
of  lower  amplitude  solitary  waves  in  both  directions  [103].  In  contrast,  a  heavy-mass 
defect  is  shifted  by  the  solitary  wave,  a  solitary  wave  is  reflected  back,  and  the 
transmitted  wave  loses  its  solitonic  character  and  is  fragmented  into  smaller  waves  of 
decreasing  amplitude  [103].  Job  and  his  collaborators  experimentally  investigated  the 
interaction  of  a  solitary  wave  with  a  light-mass  defect  in  a  ID,  unloaded  granular  crystal. 
They  showed  that  the  interaction  leads  to  the  transient  excitation  of  a  localized  mode. 
They  described  how  the  slow-timescale  local  compression  caused  by  the  solitary  wave 
around  the  defect  site  can  act  analogously  to  the  linearizing  static  compression  described 
in  the  previous  sections,  and  create  an  oscillating  localized  defect  mode  [99]. 
Starosvetsky  et  al.  also  analyzed  analytically  and  numerically  the  interaction  of  the 
solitary  wave  with  light  mass  defects.  They  used  reduced  models  that  take  into  account 
only  the  interaction  of  the  defect  mass  with  its  neighboring  particles  [104], 


6.5.4  Driven-damped  granular  crystals:  quasiperiodicity,  chaos,  and  acoustic 
rectification 

In  the  previous  sections,  we  discussed  the  existence  of  linear  and  nonlinear  localized 
modes  surrounding  defects  in  an  otherwise  homogenous  granular  crystal.  We  also 
explored  the  transient  interaction  of  traveling  solitary  waves  with  defects.  Neither  of 
these  cases  involved  a  high-amplitude  continuous  driving  force  nor  damping.  Studying 
cases  with  damping  and  continuous  driving  is  both  useful  for  real-world  applications  and 
devices,  and  involves  interesting  new  phenomena. 


In  2011,  Boechler,  et.  al.,  [58]  studied  experimentally  and  computationally  the  case  of  a 
ID  statically  compressed  granular  crystal  that  contains  a  light-mass  defect  close  to  one 


end,  and  is  subject  to  a  harmonic  driving  force,  see  left  panel  of  figure  8.  As  mentioned  in 
the  previous  section  a  light  mass  defect  will  create  a  localized  mode  with  frequency 
above  the  acoustic  band  of  the  homogenous  part  of  the  granular  crystal.  Boechler,  et.  ah, 
selected  the  frequency  of  the  driving  force  to  be  close  to  the  defect  mode  frequency. 
Because  the  driving  force  has  frequency  above  the  acoustic  band  of  the  homogenous 
granular  crystal,  the  signal  cannot  propagate  through  the  crystal  at  that  frequency. 
However,  at  sufficiently  high  amplitudes,  and  only  from  the  edge  close  to  the  defect,  a 
jump  phenomenon  occurs  from  periodic  to  quasiperiodic  and  chaotic  states,  where  the 
energy  of  the  driver  is  redistributed  to  different  frequencies  and  can  transmit  through  the 
system.  This  example  illustrates  how  the  combination  of  nonlinearity,  periodicity, 
driving,  and  asymmetric  disorder  can  create  new  material  and  device  capabilities.  In  this 
case,  this  combination  allowed  energy  to  propagate  predominantly  in  one  direction,  and 
create  an  acoustic  rectifier. 
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Figure  8:  (Left)  Schematic  diagram  of  a  ID  granular  crystal  designed  for  acoustic 
rectification  and  switching.  (Right)  Bifurcation  diagram  (Image  reproduced  from  [58], 
with  permission  from  the  publisher). 

To  understand  the  nature  of  the  bifurcations,  and  the  jump  to  the  quasiperiodic  and 
chaotic  states  that  allowed  acoustic  rectification,  Boechler,  et.  al.  conducted  parametric 
continuation  using  the  Newton-Raphson  (NR)  method  in  phase  space  [31]  and  numerical 
integration  of  the  fully  nonlinear  equations  of  motion  that  describe  the  granular  crystal. 
Dissipation  was  taken  into  account  by  using  linear  damping  (see  more  about  dissipative 
effects  in  the  next  section).  Applying  NR,  they  followed  the  periodic  family  of  solutions 
of  the  driven  system  as  a  function  of  driving  amplitude  and  studied  its  linear  stability. 
Right  panel  of  figure  8  shows  the  maximum  dynamic  force  amplitude  (four  particles  from 
the  actuator)  for  each  solution  as  a  function  of  the  driving  amplitude.  The  stable 
(unstable)  periodic  solutions  are  denoted  with  solid  blue  (dashed  black)  lines.  At  turning 
points  1,2,  stable  and  unstable  periodic  solutions  collide  and  mutually  annihilate  (saddle- 
center  bifurcation  [39]).  At  points  3,4,  the  periodic  solution  changes  stability  and  a  new 


two-frequency  stable  quasiperiodic  state  emerges  (Naimark-S  acker  bifurcation  [37]). 
Following  this  bifurcation  picture,  they  observed  in  their  experimental  setup  and 
numerical  simulations  that  with  increasing  amplitude,  a  progression  of  the  system 
response  that  followed  the  low  amplitude  stable  periodic  solution  up  to  point  1,  where  the 
system  jumps  past  the  unstable  periodic  solution  to  the  high-amplitude  stable 
quasiperiodic  state.  Further  increase  of  the  driver’s  amplitude  led  to  a  continued  cascade 
of  double  period  bifurcations  and  resulted  in  the  merging  of  distinct  frequency  peaks,  the 
fonnation  of  continuous  bands,  and  chaotic  dynamics.  As  the  quasiperiodic  and  chaotic 
states  redistribute  energy  from  the  driver  to  frequencies  within  the  transmitting  band,  it  is 
the  existence  of  these  states  which  enables  the  previously  described  acoustic  rectification. 


6.6  Dissipative  Granular  Crystals. 

Most  of  the  studies  to  date  involving  granular  crystals  ignore  dissipative  effects. 
However,  it  is  clear  from  the  experiments  that  in  many  settings  dissipation  is  strong,  and 
should  be  included.  The  sources  of  dissipation  in  granular  crystals  are  many,  including 
friction,  plasticity,  viscoelasticity,  and  viscous  drag,  among  others.  In  the  past  few  years 
there  have  been  a  number  of  analytical  and  numerical  studies  that  have  introduced 
dissipative  terms  into  the  equations  of  motion. 

In  [105],  the  authors  studied  the  effects  of  two  dissipative  mechanisms  on  the  propagation 
of  the  impulse.  The  first  was  an  intrinsic  mechanism,  an  incomplete  restitution 
mechanism  which  resulted  in  a  partial  trapping  of  the  impulse  energy  in  the  internal 
modes  of  the  grain.  The  second  mechanism  was  extrinsic,  a  velocity  dependent  friction 
/  =  -yuj .  In  both  cases,  they  showed  that  the  decay  of  the  energy  was  well  approximated 
by  an  exponential  function.  The  attenuation  of  traveling  pulses  in  ID  unloaded  granular 
crystals  due  to  on-site  linear  damping  /  =  —yut ,  was  also  analyzed  in  [106].  They  found 

an  overall  exponential  decay  of  the  energy,  which  depends  on  the  exponent  of  the 
interaction  potential,  and  causes  the  pulse  to  slow  down  as  it  moves.  They  also  showed 
that  the  shape  and  the  width  of  the  pulse  remain  unchanged. 

Job  and  his  collaborators  studied  the  interaction  of  a  solitary  wave  with  boundaries  in  a 
ID  granular  crystal,  considering  two  dissipative  mechanisms:  internal  viscoelasticity  and 
solid  friction  of  the  beads  due  to  their  weight  on  the  track.  Viscoelastic  dissipation  was 
taken  account  by  considering  a  dissipative  force  at  the  contact  of  the  two  beads  in  the 
fonn  /  =  i]Adl([u:  ]  -  «,]f2  )  [  1 07],  where  /;  includes  unknown  coefficients  due  to 
internal  friction  of  the  material.  Solid  friction  was  included  by  considering  a  force 
/  =  jumg .  As  a  result  of  these  terms,  the  dissipation  is  expected  to  produce  broader 
solitary  waves. 


In  [108]  viscous  dissipation,  depending  on  the  relative  velocity  between  neighboring 

particles  was  included  in  the  model  as  /  =  p(«M-2M;+«i+1),  where  p  is  the  viscosity 
coefficient.  The  authors  investigated  its  influence  on  the  shape  of  a  steady  shock  wave.  In 
[109],  they  solved  the  following  system  of  nondimensional  equations: 

ii,  =  [p(ui+ 1  -«,■)-  0,  ~  um  )"  W(ui  ~  uM )  -  u, )  +  (u,  ,  -  Uj )"  ]0(uM 

where  9  is  the  Heaviside  function.  They  found  the  inclusion  of  this  relative  velocities 
dependent  viscous  damping  may  yield  interesting  effects  such  as  the  creation  of 
secondary  pulses.  A  different  approach  was  presented  in  [110],  where  the  authors 
provided  a  quantitative  characterization  of  dissipative  effects  for  the  solitary  wave 
propagation  in  ID  granular  crystal.  They  incorporated  a  phenomenological  nonlinear 
dissipation  that  depends  on  the  particle’s  relative  velocities.  By  using  optimization 
schemes  and  experiments,  they  calculated  a  common  dissipation  exponent  with  a 
material-dependent  prefactor. 

Most  of  the  above  studies  concern  the  attenuation  of  propagating  pulses  generated  by  an 
impulsive  excitation.  Recent  experiments  in  ID  compressed  granular  crystal,  subject  to  a 
continuous  harmonic  driving  at  one  end,  also  revealed  a  strong  attenuation  of  the  signal 
[58].  To  account  for  the  dissipation  in  these  experimental  settings,  we  used  a  linear  on¬ 
site  damping  /  =  —yui  with  a  damping  coefficient  y  that  was  selected  to  match  the 
experimental  results. 


6.7  Two-dimensional  granular  crystals. 

The  richness  of  the  nonlinear  dynamic  phenomena  found  in  one-dimensional  systems 
points  at  the  study  of  higher  dimensional  systems  for  the  discovery  of  new  dynamic 
effects.  For  example,  two-  and  three-dimensional  nonlinear  systems  are  expected  to 
present  additional  families  of  wave  modes  not  realizable  in  the  ID  case;  new  types  of 
solitary  waves  propagating  in  the  axial  and  lateral  directions  (particularly  interesting  for 
wave  energy  redirection  and  wave  guiding);  complex  nonlinear  resonance  interactions 
occurring  between  spatially  extended  modes  and  localized  waves;  and  enhanced 
possibilities  for  acoustic  wave  energy  localization  and  trapping  across  spatial  or  temporal 
scales. 

The  dynamic  properties  of  2-D  granular  crystals  have  only  been  partially  characterized. 
In  particular,  experimental  efforts  are  few  and  sparse  although  such  systems  are  expected 
to  present  a  rich  variety  of  novel  dynamic  phenomena.  In  the  literature,  several  authors 
proposed  different  models  to  characterize  the  mechanical  response  of  two-dimensional, 
ordered  granular  media.  For  example,  [111]  presented  a  model  for  a  square  lattice  of 
elastically  interacting  particles,  including  relative  particles’  rotation.  Toumat  et  ah,  2011 


[112]  proposed  a  theoretical  model  to  describe  out-of-plane  elastic  waves  in  a  monolayer 
granular  membrane  consisting  of  a  hexagonal  lattice  of  particles.  Their  model  was  the 
first  one  to  include  shear  and  bending  rigidity  at  the  contact  between  particles,  and  to 
calculate  dispersion  relations  reporting  the  presence  of  acoustic  band  gaps. 

The  simplest  example  of  a  highly  nonlinear  2-D  granular  crystal  consists  of  a  unifonn, 
uncompressed  square  packing  of  elastic  particles  in  contact  with  each  other.  When  this 
system  is  excited  on  one  side  by  a  unifonn,  planar  waveform,  its  response  is  expected  to 
be  quasi-one-dimensional  [13]  and  the  response  of  the  system  can  be  characterized  by  a 
“curtain  solution”  derived  similarly  to  equation  (10)  in  this  chapter.  The  first 
experimental  characterization  of  the  dynamic  behavior  of  a  square  packing  of  particles 
was  provided  in  [81],  using  photoelastic,  elliptical  disks,  excited  by  an  explosive  charge. 
The  same  study  characterized  the  stress  propagation  response  of  various  geometrical 
packing  of  elliptical  disks,  and  concluded  that  the  contact  nonnals  and  the  vectors 
connecting  particles'  centers  of  mass  influence  the  wave  propagation  characteristics,  such 
as  load  transfer  path  and  load  attenuation.  Discrete  element  numerical  models  (DEM) 
were  used  to  analyze  similar  systems  [113]. 

The  presence  of  solitary  waves  forming  and  traveling  in  2-D  square  granular  crystals  was 
reported  and  studied  quantitatively  for  the  first  time  by  [114]  using  triaxial 
accelerometers  embedded  within  selected  particles  in  the  crystals.  A  larger  number  of 
studies  focused  on  the  dynamic  behavior  of  hexagonal  packing  under  different  loading 
conditions  [81,  113,  115-121],  One  of  the  major  difficulties  in  the  experimental 
realization  of  acoustic  materials  based  on  two-dimensional  nonlinear  granular  lattices  is 
the  sensitivity  of  such  systems  to  the  presence  of  inaccuracy  of  the  particles  geometry.  In 
the  ideal  configuration,  all  particles  have  an  equal  number  of  contacts  and  equal 
equilibrium  forces.  The  presence  of  small  defects  in  experiments,  however,  can  lead  to 
the  loss  of  contact  between  particles  or  to  the  local  compression  in  the  surrounding 
particles.  Such  loss  of  contacts  or  local  compression  ultimately  results  in  a  disordered 
energy  transfer  between  the  particles.  A  few  past  works  studied  the  effects  of 
imperfections  in  two-dimensional  granular  crystals  and  their  role  in  the  stress  wave 
propagation  [115-118].  While  Herztian  behavior  predicts  a  1/6  power-law  between 
maximum  force  and  wave  speed  [46],  it  was  found  that  the  presence  of  defects  tends  to 
increase  the  wave  propagation  speed  to  a  1/4  power  law  relationship,  effectively  inducing 
deviations  from  the  theoretical  Hertzian  behavior.  This  deviation  from  Hertzian  behavior 
was  observed  only  for  granular  crystals  with  low  precompression.  Increasing  the 
precompression  applied  on  hexagonal  arrays  was  seen  to  cause  a  transition  to  a  fully 
Hertzian  behavior  [115-118],  More  recently  Leonard  et  al.  characterized  experimentally 
the  response  of  regular  2D  square  granular  crystals  and  showed  that  variation  in  the 
packing  geometry/composition  (figure  9)  can  dramatically  vary  the  directional  wave 
propagation  in  such  systems  [122], 
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Figure  9:  (Left)  Dynamic  response  of  a  two-dimensional  granular  crystals  formed  by  a 
square  centered  packing  of  cylinders  and  spheres  of  different  materials  (see  inset). 
Variation  of  the  materials  configuration  lead  to  dramatic  changes  of  the  wave  propagation 
front,  as  shown  from  experiments  and  numerical  simulations  [122].  (Right)  Design 
concept  of  a  tunable,  nonlinear  acoustic  lens  obtained  with  a  two-dimensional  array  of 
particle  chains.  The  formation  of  the  focal  spot  (i.e.,  the  “sound  bullet”)  is  evident  on  the 
host  medium  on  the  right. 

Two-dimensional  arrays  of  particles  have  also  been  shown  to  fonn  tunable  acoustic 
lenses  (figure  9,  right),  supporting  the  formation  of  concentrated  acoustic  pulses  at  the 
focal  point  (“sound  bullets”,  [57]).  The  fundamental  ability  to  redirect  nonlinear  acoustic 
pulses  in  two-dimensional  systems  has  also  been  studied  looking  at  the  pulse  splitting  and 
recombination  in  a  y-shaped  granular  network  [123-125].  These  works  showed 
theoretically,  numerically  and  experimentally  the  ability  to  bend  and  split  incident  pulses, 
and  redirect  mechanical  energy  as  a  function  of  the  branch  geometry. 

Additional  work  on  the  dynamic  behavior  of  ordered  two-dimensional  granular  crystals  is 
needed  to  fully  understand  the  dynamic  response  of  such  systems,  and  to  characterize 
how  these  properties  depend  on  the  underlying  particle  arrangement.  Variations  of  the 
excitation  type  (impulsive  or  harmonic  forcing)  is  expected  to  lead  to  the  discovery  of 
interesting  new  acoustic/dynamic  phenomena  including  wave  guiding,  trapping,  filtering 
and  localized  breathing  modes. 

6.8  Future  Directions  and  Conclusions 

The  preceding  chapters  of  this  book  have  demonstrated  how  structural  periodicity  can  be 
utilized  to  create  new  materials  with  unprecedented  physical  properties.  In  such  materials 
the  individual  building  blocks  are  assembled  in  carefully  designed  structures,  where  by 
working  together,  they  cause  the  bulk  material  to  present  properties  greater  than  those  of 
the  individual  components.  This  general  concept  of  obtaining  “materials  by  design”  is  not 
new,  and  has  been  a  long-term  quest  for  chemists  and  material  scientists  alike.  For 


instance,  chemists  have  long  been  trying  to  engineer  crystals  and  molecules  by  arranging 
atoms  in  specific  lattices  and  geometries,  to  obtain  a  specific  bulk  property.  However,  by 
extending  this  concept  past  molecules  and  crystal  grains,  to  specially  designed  structural 
building  blocks  -  from  the  nano  to  macroscales  -  we  open  a  whole  new  field  of 
possibilities. 

One  of  the  main  benefits  of  such  designed  materials  is  that  they  enable  new  technological 
capabilities.  In  terms  of  structural  materials,  new  materials  with  multifunctional 
properties  become  available,  where  they  have  both  structural  and  additional  dynamic 
properties.  Perhaps  more  importantly,  by  creating  materials  with  new,  previously  unseen 
properties,  new  devices  are  enabled.  By  using  these  designed  materials  in  devices, 
engineers  have  new  capabilities,  and  can  create  unprecedented  applications.  Furthermore, 
as  such  materials  are  “designed”  by  construction,  they  can  be  tailored  for  use  in 
specifically  targeted  applications. 

The  range  of  possible  bulk  responses  from  such  designed  materials  depends  in  part  on  the 
complexity  of  the  interaction  of  the  fundamental  building  blocks  -  such  as  the  particles  in 
granular  crystals  or  resonators  in  metamaterials.  As  described,  the  design  of  these 
periodic  structures  was  initially  based  on  linear  interactions.  The  presence  of  nonlinearity 
in  these  systems  is  an  added  advantage.  In  this  chapter,  we  predominately  focused  on 
nonlinear  dynamic  phenomena  in  granular  crystal  systems,  where  the  nonlinearity  was 
caused  by  the  geometric  inter-particle  interactions  between  elastic  particles.  As  described, 
nonlinear  dynamics  causes  the  existence  of  new  useful  dynamic  phenomena  and  coherent 
structures.  This  includes  solitary  waves,  discrete  breathers,  bifurcations,  quasiperiodicity, 
and  chaos,  among  others.  Nonlinearity  enables  a  dramatic  tunability  of  the  material 
responses,  by  providing  an  unprecedented  sensitivity  to  variations  of  materials  and 
external  parameters.  This  tunability  also  creates  new  ability  to  scale  the  system’s 
components  and  overall  size,  by  breaking  the  frequency-discreteness  restrictions. 

Because  of  the  inherent  complexity  of  nonlinear  systems,  which  paradoxically  enables 
such  useful  phenomena,  the  richness  of  the  dynamic  phenomena  expected  by  nonlinear 
periodic  acoustic  systems  and  their  potential  applicability  to  new  engineering  systems  is 
still  mostly  unexplored.  The  ability  to  obtain  and  understand  such  systems  will  rely 
heavily  on  the  ability  to  develop  predictive  theoretical  and  computational  tools,  to  guide 
the  design  and  testing. 

Some  future  areas  of  interest,  with  respect  to  the  study  of  nonlinear  periodic  phononic 
structures,  include,  but  are  not  limited  to  the  following.  As  nonlinearity  has  been  applied 
to  spring-mass-like  systems,  in  granular  crystals,  nonlinearity  could  further  be  applied  to 
the  study  of  nonlinear  metamaterials,  nonlinear  resonant  structures,  or  superlattices  with 
nonlinear  elastic  components.  The  study  of  hybrid  linear-nonlinear  systems,  in  which  the 
interaction  of  linear  modes  and  deformations  with  nonlinear  effects  could  lead  to  the 


observation  of  new  dynamic  phenomena  like,  for  example,  the  amplitude  dependent 
filtering  of  acoustic  signals  [126]. 

The  design  of  new  material  systems  where  nonlinear  material  responses  interplay  with 
active  building  blocks  or  other  multi-physical  effects,  like,  for  example,  the  use  of 
ferroelectricity  to  tune  the  stiffness  of  materials,  is  yet  another  area  that  could  lead  to  the 
discovery  of  unprecedented  material  response  relying  on  the  conversion  of  mechanical  to 
electrical  energy.  The  ability  to  couple  multi-physical  effects  in  materials  with  the 
creation  of  periodic  structures  can  also  lead  to  the  creation  of  tunable  multifunctional 
devices,  like  novel  opto-mechanical  sensors  [127],  or  phoXonic  systems  [128],  For 
example,  the  generation  of  nonlinear  modes  in  nonlinear  acoustic  crystals  could  lead  to 
accurate  frequency  conversion,  and  the  presence  of  breathing  modes  or  nonlinear 
localized  modes  can  be  exploited  for  energy  localization  and  harvesting. 

Because  of  the  similarity  of  acoustic  and  elastic  wave  propagation  to  phonon 
propagation,  through  scaling,  the  effects  studied  here  could  also  be  extended  to  smaller 
scales,  which  affect  the  propagation  of  heat.  For  instance,  as  described  in  acoustics, 
nonlinear  periodic  structures  have  been  utilized  to  create  tunable  rectifiers  based  on  the 
onset  of  bifurcation  instabilities,  which  could  provide  new  ways  to  control  the  flow  of 
acoustic  energy,  create  acoustic  logic  devices,  and  be  used  in  novel  energy  harvesting 
systems  [58],  These  same  ideas  could  be  scaled  down  to  create  new  ways  to  control  heat 
propagation,  create  materials  with  direction-dependent  thennal  conductivities,  and  create 
thermal  logic  devices. 

Furthermore,  the  newly  explored  capabilities  of  granular  crystals,  and  nonlinear  periodic 
phononic  structures,  should  be  further  explored  for  their  potential  in  engineering 
applications.  The  ability  to  engineer  the  dispersion  relation  through  nonlinearities  could 
be  implemented  in  tunable  vibration  filtering  devices  and  in  noise  and  vibration 
insulating  systems.  Compact  solitary  waves  with  robust  properties  and  large  amplitudes 
could  find  use  in  biomedical  devices  with  improved  resolution  and  signal-to-noise  ratio 
[57],  or  in  the  non-destructive  evaluation  of  materials  [56]. 

The  study  of  nonlinearity  in  engineered  materials,  like  phononic  crystals  and 
metamaterials  is  still  at  an  early  stage  of  development.  By  understanding  the  fundamental 
properties  of  nonlinear  acoustic  crystals,  nonlinear  phononic  systems  and  nonlinear 
resonant  structures,  new  physical  phenomena  can  be  discovered,  and  lead  to  a  new  ability 
to  design  and  implement  materials  and  devices. 
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